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WAVELENGTH DIVISION MULTIPLEXED 
OPTICAL SOLITONS 

CROSS-REFERENCE TO RELATED 
APPLICATIONS 

This application claims benefit of U.S. Provisional appli- 
cation No. 60/153,276, filed Sep. 9, 1999. 
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The fundamental equations governing M numbers of 
co-propagating waves in a nonlinear fiber including the 
CPM phenomenon are the coupled nonlinear Schrodinger 
equations 
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STATEMENT AS TO FEDERALLY-SPONSORED 
RESEARCH 

The invention described herein was made in the perfor- 
mance of work under a NASA contract, and is subject to the 
provisions of Public Law 96-517 (U.S. C. 202) in which the 
contractor has elected to retain title. 


(7=1,2, 3,... M) 
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Here, for the jth wave, Ay(z,t) is the slowly-varying 
amplitude of the wave, v -, the group velocity, (3 2j -, the 
dispersion coefficient (P 2 /=dv^- _1 /da)), a p the absorption 
15 coefficient, and 


BACKGROUND 

A soliton is a light pulse that maintains its shape over a 
long distance. A soliton wave does not broaden, does not 
loose its shape and does not substantially weaken as it 
travels across a particular medium. This means that better 
transmission and fewer repeaters are required. Solitons have 
been known to carry 5 gigabits of information per second 
over 15,000 kilometers or 10 gigabits per second over 
11,000 kilometers. 

Different ways have been known to create the soliton in 
a fiber. 


.. = n ^J (2) 

J cAeff 

is the nonlinear index coefficient with A eff as the effective 
core area and n 2 =3.2xl0 -16 cm 2 /W for silica fibers, w } - is the 
carrier frequency of the jth wave, c is the speed of light, and 
z is the direction of propagation along the fiber. 
Introducing the normalizing coefficients 
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However, wavelength division multiplexed systems have 
provided a challenge for soliton communications. In wave- 30 
length division multiplexing, multiple beams are provided. 
Each of the beams has its own wavelength. This produces, 
however, complex non-linear interaction between the propa- 
gating pulses on these beams. This non-linear interaction has 
hampered the existence of solitons in WDM beams in fibers. 35 

Attempts have been made to compensate for the effects. 

SUMMARY 

The present application teaches modification of various 40 
parameters of an optical system in order to produce solitons 
for optical communications through an optical fiber. The 
operation is done by setting the parameters of the pulses on 
these beams that will be launched into the fibers according 
to a specified relationship. 45 

BRIEF DESCRIPTION OF THE DRAWINGS 

These and other aspects will now be described in detail 
with reference to the accompanying drawings, wherein: 

FIG. 1 shows a fiber with a plurality of WDM solitons 50 
therein. 

DETAILED DESCRIPTION 

Solitons are based on a delicate balance between disper- 55 
sion effect and the non-linear effect. A temporal soliton is a 
specially shaped optical pulse which can propagate in a fiber 
without shape degradation. 

When beams with different wavelengths are 
co-propagating in a single mode fiber 100 as shown in FIG. 60 
1, such as in a wavelength division multiplexed fiber system, 
interaction among the pulses of the different beams can 
cause non-linear cross phase modulation. This can effect the 
integrity of solitons in wavelength division multiplexed 
beams such at 110 and 115. Different techniques have been 65 
attempted to prevent this; often called “dispersion manage- 
ment”. 
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a =1,2, 3,... M) 


Here, T 0 is the pulse width, P 0j - is the incident optical 
power of the jth beam, and d l7 -, the walk-off parameter 
between beam 1 and beam j, describes how fast a given pulse 
in beam j passes through the pulse in beam 1. In other words, 
the walk-off length is 

Lw(ij) = ^of\dij\- ( 7 ) 

So, L wQjy is the distance for which a faster moving pulse 
(say, in beam j) completely walked through the slower 
moving pulse in beam 1. The nonlinear interaction between 
these two optical pulses ceases to occur after a distance 
L W ( l7 ->. For cross phase modulation to take effect 
significantly, the group -velocity mismatch must be held to 
near zero. 

Finding the analytic solution of Eq. (6), which is a set of 
simultaneous coupled nonlinear Schrodinger equations, may 
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be a formidable task. However, it may be solved numerically 
by the split-step Fourier method, which was used success- 
fully earlier to solve the problem of beam propagation in 
complex fiber structures, such as the fiber couplers, and to 
solve the thermal blooming problem for high energy-laser 5 
beams. According to this method, the solutions may be 
advanced first using only the nonlinear part of the equations, 
and then the solutions are allowed to advance using only the 
linear part of Eq. (6). This forward stepping process is 
repeated over and over again until the desired destination is 10 
reached. The Fourier transform is accomplished numerically 
via the well-known fast Fourier transform technique. 

It is well know that, for an idealized fiber with no loss, an 
optical soliton on a single wavelength beam takes the initial 
form 15 

u(0,x)=N sech(t) (8) 

where N is the soliton magnitude and 

n 2 _ ( 9 ) 20 
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It is also known that the single beam soliton equation is 
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it may be possible to construct a set of initial pulses which 
will propagate in the same manner as the single soliton pulse 
case, i.e., undistorted and without change in shape for 
arbitrarily long distances. If the initial pulses are selected as 
follows: 

Uj(0 ,x) = (1 +2(M- sech(x) 

0=1, 2, 3, . . . M), (12) 

where M is the number of WDM beams. 

Using these initial pulse forms numerical simulation is 
carried out to solve Eq. (11). The split-step Fourier method 
is used. The fiber parameters used for the situation are L, 
length of fiber equal to 1000 km; p 2 =dispersion coefficient, 
equal to -2 ps 2 /km; y is the nonlinear index coefficient, equal 
to 20 W -1 km -1 ; T 0 , pulse width, equal to 10 ps; L z> =50 km; 
and L z?1 =50 km. 

Four cases with M=l, 2, 3, 4 are treated. The M=1 case 
corresponds to the well-known single soliton case; here, the 
amplitude for the fundamental soliton is 1. For the two-beam 
case, the amplitude is (3) _y2 =0.57735. For the three-beam 
case, it is (5) -y2 =0.4472136. For the four-beam case, it is 
(7) _y2 =0.37796447. In other words, using Eq. (5), the initial 
power for M beams can be found: 

m i (i3) 

0j r T§[l+2(M-l)] 


Here, the dispersion length L D and the nonlinear length 
are defined earlier in Eq. (6). In the case of anomalous 30 
group velocity dispersion (GVD) for soliton, sgn((3 2 )=-l. 

For the fundamental soliton case, N=l. This means that 
when an initial pulse with pulse shape given by Eq.(8), with 
a normalized amplitude of unity, is launched inside an ideal 
lossless fiber, the pulse will retain its hyperbolic secant 35 
shape without degradation for arbitrarily long distances, e.g., 
over thousands of kilometers. The delicate balance between 
the dispersion effect represented by L D and the nonlinear 
self phase modulation effect represented by occurs at 
N=1 for the fundamental soliton. The nonlinear effect on a 40 
pulse for a single wavelength beam is embodied in L nl , 
while the dispersion effect on the pulse is embodied in L D . 

An embodiment considers the case of multiple wave- 
length division multiplexed (“WDM”) beams in a fiber. 
Starting with an idealized fiber which is lossless (i.e., 0^=0 45 
for all beams) and which possesses uniform group velocity 
dispersion (i.e., v s -=v^ for all beams) within the wavelength 
range under investigation, the equations governing the 
propagation characteristics of signal pulses are: 
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(y = l,2, 3,... M). 
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The anomalous GVD case in which sgn(p 2 )=-l is con- 
sidered. It is seen from the above equation that the summa- 
tion term representing the cross phase modulation effect is 
twice as effective as the self-phase modulation (SPM) effect 60 
for the same intensity. This observation also provides the 
idea that cross phase modulation may be used in conjunction 
with self-phase modulation on the WDM pulses to counter- 
act the GVD effect, thus producing WDM solitons. Com- 
paring the bracketed terms in Eqs. (10) and (11) shows that 65 
if one chooses the correct amplitudes for the initial pulses on 
WDM beams and retains the hyperbolic secant pulse form, 


The amplitude of the fundamental solitons on WDM 
multibeams becomes successively smaller as the number of 
beams is increased. This is because the nonlinear effect 
becomes more pronounced when more beams are present. 
Numerical simulation shows that after propagating 1000 km 
through this fiber the original pulse shape for all these WDM 
pulses remains unchanged. 

As additional examples, if a silica fiber is used at 1.55 um, 
with a 1 ps pulse, a P of 5 watts is obtained. With the same 
materials, but a 10 ps pulse, a 50 mw pulse is obtained. 

When this is used, normalized power is said to be “1”. The 
anomalous group velocity dispersion case in which sgn 
((3 2 y)=— 1, is considered. It is seen from the above equation 
that the summation term in the bracket representing the cross 
phase modulation (CPM) effect is twice as effective as the 
self phase modulation (SPM) effect for the same intensity. 
Hence, cross phase modulation may be used in conjunction 
with self phase modulation on the WDM pulses to counter- 
act the group velocity dispersion effect. WDM solitons are 
thus produced. 

Comparing the bracketed terms in Equations (10) and (11) 
shows that if one chooses the correct amplitudes for the 
initial pulses on WDM beams and retains the hyperbolic 
secant pulse form, it may be possible to construct a set of 
initial pulses which will propagate in the same manner as the 
single soliton pulse case, i.e., undistorted and without 
change in shape for arbitrarily long distances. Let us choose 
the initial pulses as follows: 

Wy(0,T)=(l+2(M-l) -% sec h(z t) 

0=1, 2, 3, ... M), (12) 

where M is the number of WDM beams. 

Using these initial pulse forms numerical simulation was 
carried out to solve Equation (11). The split-step Fourier 
method was used. The fiber parameters used for the simu- 
lation are: 

L=length of fiber=1000 km 
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p 2 =clispersion coefficient =-2 ps 2 /km 
y=nonlinear index coefficient=20 w _1 km _1 
T 0 =pulse width=10 ps. 

1^=50 km 5 

km. Four cases with M=l, 2, 3, 4 were treated. 
The M=1 case corresponds to the well known single 
soliton case; here, the amplitude for the fundamental 
soliton is 1. For the 2-beam case, the amplitude is 
(3) _V2 =0. 57735. For the 3-beam case, it is (5)“ V2 = 10 

0.4472136. For the 4-beam case, it is (7) -y2 = 
0.37796447. It is noted that the amplitude of the 
fundamental solitons on WDM multi-beams becomes 
successively smaller as the number of beams is 
increased. This is because the nonlinear effect becomes ^ 
more pronounced when more beams are present. 

Another way of looking at the soliton is that. the correct 
power to form a soliton is obtained by the value P; 

p_ & 20 

ygTl 

where beta 2 is the dispersion coefficient, which is a number 
that depends on the fiber; gamma is the non linear coefficient 
of the fiber material, and T 0 is proportional to the pulse 25 
width. 

Numerical simulation shows that after propagating 1000 
km through this fiber the original pulse shape for all these 
WDM pulses remains unchanged. It thus appears that the 
initial forms chosen for the pulses on WDM beams are the 30 
correct soliton forms for WDM beams. 

What is claimed is: 

1. A method, comprising: 

forming a plurality of wavelength division multiplexed 35 
hyperbolic secant shaped optical beams; and 
simultaneously modulating all of said plurality of secant 
shaped optical beams as wavelength division multi- 
plexed beams in a hyperbolic secant form, and all of 
said beams being initially modulated with a specified 40 
effective power that forms wavelength division multi- 
plexed solitons. 

2. A method as in claim 1 wherein said modulating 

comprises determining a single effective power for said] 
optical beams which allows said solitons to exist. 45 

3. A method, comprising 

forming a plurality of wavelength division multiplexed 
hyperbolic secant shaped optical beams; and 
modulating said secant shaped optical beams as wave- 
length division multiplexed beams in a hyperbolic 50 
secant form, and said beams being initially modulated 
with a specified effective power that forms wavelength 
division multiplexed solitons, wherein said modulation 
comprises determining a single effective power for said 
optical beams which allows said solitons to exist, 55 
wherein said amplitude is selected such that a self 
phase modulation effect and cross-phase modulation 
effect counteract a group velocity dispersion effect to 
form said solitons which will propagate in a fiber for at 
least 10000 km. 60 

4. A method as in claim 2, wherein there are M wave- 
length division multiplexed beams, and said effective power 
is a function of M where M is a number of said multiplexed 
beams, M being an integer greater than 1. 

5. A method as in claim 4, wherein there are M=2 beams, 65 
and the amplitude is 1/^3 times the square root of the 
specified power level for a single beam soliton. 
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6. A method as in claim 4, wherein the specified amplitude 
level is 

1 

Vl +2(M-1) 

times the square root of the specific power level for a single 
beam soliton, where M is the number of beams, and M>1. 

7. A method of conducting beams, comprising: 
choosing a series of initial pulses, each of which have a 

hyperbolic secant shape, and which have effective 
powers which are set to allow soliton transmission of 
all of the serves of pulses, in a wavelength division 
multiplexed form, wherein the specified amplitude 
level is 

1 

Vl +2(M-1) 

times the square root of the specific power level for a 
single beam soliton, where M is the number of beams, 
and M>1. 

8. A method of keeping parallel beams of information 
from degradation when propagating through an optical fiber 
comprising: forming WDM solitons which has a specified 
amplitude level is 

1 

Vl +2(M-1) 

times the square root of the specific power level P for a 
single beam soliton, where M is the number of beams, and 
M>1, and 



which pulses have approximately a hyperbolic secant ampli- 
tude envelope and which also have a pulse width x and a 
peak power related where p 2 is the dispersion coefficient, 8 
is the nonlinear coefficient of the fiber material, and T 0 is a 
value that is proportional to the pulse width. 

9. A method, comprising 

forming a plurality M of wavelength division multiplexed 
beams on a single fiber, which have powers that are set 
to have uniform group velocity dispersion, and which 
are related such that the sign of (p 2; -)=-l, and which 
have an initial form u / -(0,t)=^ 1 +2(M-l)sec h(x), where 
j is an integer<=M, and M is an integer >1, (3 is the 
dispersion coefficient, and t is the pulse width. 

10. A method comprising: 

forming a plurality of wavelength division multiplexed 
shaped co-prop agating parallel pulses on wavelength 
division multiplexed (WDM) optical beams; and 

modulating said optical beams as wavelength division 
multiplexed beams co -prop agating pulses in a soliton 
form, and said wavelength division multiplexed beams 
being initially modulated with a specified effective 
power that forms wavelength division multiplexed soli- 
tons and which power is dependent on a number of said 
pulses. 

11. A method as in claim 1, wherein said modulating 
comprises determining the effective power for said all of 
said beams, which allows said WDM solitons to exist. 



US 6,612,743 B1 


7 

12. A method as in claim 10, wherein said amplitude for 
the co-propagating bit-parallel pulses is selected such that a 
self phase modulation effect and cross-phase modulation 
effect counteract a group velocity dispersion effect to form 
said solitons which will propagate in a fiber for at least 5 
10000 km. 

13. A method as in claim 11, wherein there are M 
co-propagating beams which form bit-parallel pulses on M 
wavelength division multiplexed beams, and said effective 
power is a function of M, where M is an integer >1 
representing the number of co-propagating bit-parallel 
pulses. 

14. A method as in claim 13, wherein there are M=2 
co-propagating bit-parallel pulses on 2 WDM beams, and 
the amplitude is 1/^3 times the square root of the specified 
power level for a single beam soliton. 

15. A method as in claim 13, wherein the specified 
amplitude lever is 

1 

Vl +2(M-1) 

times the square root of the specific power level for a single 
beam soliton, where M is the number of co-propagating 
bit -parallel pulses on WDM beams, and M>1. 

16. A communication system, operating to transmit bit- 
parallel data, comprising: 

an optical cable, and 

an optical modulating element forming a plurality of 
wavelength division multiplexed hyperbolic secant 
shaped co-propagating pulses on wavelength division 
multiplexed (WDM) optical beam and modulating said 
secant shaped WDM optical beams as wavelength 
division multiplexed beams in a hyperbolic secant form 
with a specified effective power that forms wavelength 
division multiplexed solitons, 

wherein said modulating element selects an amplitude for 
the co-propagating parallel pulses is selected such that 
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a self phase modulation effect and cross-phase modu- 
lation effect counteract a group velocity dispersion 
effect to form said solitons which will propagate in a 
fiber for at least 10000 km. 

17. A communication system, operating to transmit bit- 
parallel data, comprising: 

an optical cable, and 

an optical modulating element forming a plurality of 
wavelength division multiplexed hyperbolic secant 
shaped co-propagating pulses on wavelength division 
multiplexed (WDM) optical beam and modulating said 
secant shaped WDM optical beams as wavelength 
division multiplexed beams in a hyperbolic secant form 
15 with a specified effective power that forms wavelength 
division multiplexed solitons, 
wherein said modulating element producing M 
co-propagating bit-parallel pulses on M wavelength 
division multiplexed beams, and said effective power is 
a function of M, where M is the number of 
co-propagating bit-parallel pulses and M>1. 

18. A system as in claim 17, wherein said modulating 
element producing M=2 co-propagating parallel pulses on 2 
WDM beams, and the amplitude is 1/^3 times the square 
root of the specified power level for a single beam soliton. 

19. A system as in claim 18, wherein the specified 
amplitude level is 

1 

Vl +2(M-1) 

times the square root of the specific power level for a single 
beam soliton, where M is the number of co-propagating 
35 bit-parallel pulses on WDM beams, and M>1. 





